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Abstract 



jis ■ We present formulae for the mass differences AM^ and AMg in the B^ s'^d s ^y^' 

D . tems and for the CP violation parameter e which are valid in minimal flavour viola- 

tion models giving rise to new four-fermion AF = 2 operators. Short distance con- 
tributions to AMg , AMrf and e are parameterized by three real functions Ff-/. , F^ and 
X ; Fit, respectively (F/^ = F^ = F^^ holds only if the Standard Model {V-A)(g){V- A) 

C^ ' operators dominate). We present simple strategies involving the ratio AMg/AMd, 

sin 2/3 and 7 that allow to search for the effects of the new operators. We point out 
that their sizable contributions to the ratio AMg/AMa would in principle allow 7 
to be larger than 90°. Constraints on the functions Fj\ imposed by the present (and 
future) experimental data are also discussed. As an example we show that for large 
tan/3 = v^jvx and iJ"*" not too heavy, Fj^ in the MSSM with heavy sparticles can 
be substantially smaller than in the SM due the charged Higgs box contributions 
and in particular due to the growing like tan^ ^ contribution of the double penguin 
diagrams involving neutral Higgs boson exchanges. As a result the bounds on the 
function F^^ can be violated which allows to exclude large mixing of stops. In this 
scenario the range of sin 2/3 following from e and AM^ is identical to the SM ones 
(0.5 < sin 2/3 < 0.8). On the other hand 7 following from AM^j AM^ is lower. 



1 Introduction 

The determination of the CKM parameters and of the related unitarity triangle (UT) 
are the hot topics in particle physics [|l|. In this context a clean measurement of the 
angle j3 in the unitarity triangle through the time dependent CP asymmetry, a^i^g(t), in 
B^{B^) — >• ipKs decays is very important. 
In the Standard Model (SM) 

a^Ks{t) = - a^Ks sin(AMrft) = - sin2/3sin(AMdt), (1.1) 

thereby allowing direct extraction of sin 2(3. The most recent measurements of a^Xg from 
the BaBar and Belle Collaborations give 

_ J 0.59 ±0.14 ±0.05 (BaBar) [2] 
"^^' ~ I 0.99 ±0.14 ±0.06 (Belle) §. ^ ' ^ 

Combining these results with earlier measurements at CDF (0.79tQ44) and by the 
ALEPH collaboration (0.84t5;o4 i 0.16) p would give the grand average 

a^^3 = 0.79 ±0.10 , (1.3) 

but in view of the fact that BaBar and Belle results are not fully consistent with each 
other we believe that a better description of the present situation is a^Ks = 0-80 i 0.20. 
Similarly important for the determination of the unitarity triangle will be the mea- 
surement of the ratio of the mass differences in the B^g-B^^ systems, AMg/AMd. The 
experimental values of AM^ and AM^ read P] 

AMd = (0.487 ± 0.009)/ps , AM, > 15.0/ps , (1.4) 

implying 

AM, , , , 

which is compatible with the SM expectations. Because theoretically this ratio is con- 
siderably cleaner than AM, and AM^ themselves, its precise measurement will have an 
important impact on the determination of the unitarity triangle and on the tests of the 
SM and its various extentions. As emphasized in [^ |^, H, ^ [10| this impact will be 
even stronger in conjunction with the measurement of a^Xs- It is therefore exciting that 
AMg/AMd should be measured already this year in Run II at Fermilab, while improve- 
ments of the a^Ks measurements are expected from BaBar, Belle, CDF, DO and, at later 
stages, from BTeV and LHC experiments. 
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The result (p.. 3D for a^Ks should be compared with the value of sin 2(3 obtained from the 



analyses of the unitarity triangle in the framework of the Standard Model (SM) [jl], |7|, |ri[] 
that center around (sin2/5)sM ~ 0.70 with estimated errors ranging from 0.07 to 0.24. 
Clearly in view of large experimental error in ( p. .31) and the considerable uncertainty in 
the error estimates of (sin 2/5) sm, the SM fits are compatible with the experimental value of 
a^Ks ■ The small value of a^xs found earlier by the BaBar collaboration |jl2| has not been 



confirmed by most recent data. On the other hand, the large value of a^xs measured by 
the Belle collaboration may be suggestive of new physics contributions to B^{B^) —>■ ipKs, 
B^g-B^g mixing, K^-K^ mixing and/or the parameter e measuring CP violation in the 
K^-K^ mixing. Thus the analyses [jTU], |]T^-|^ of new contributions done in the context 



of small a^Ks value reported by BaBar |T2[ could still be relevant if properly reformulated. 
Such new contributions could modify not only the relation between a^Ks ^i-nd sin 2/5 in 
( |1 . ID but also the value of sin 2/5 obtained from the fits of the unitarity triangle. 

In general models of new physics potentially contributing to a^Ks i ^^d,s and/or e fall 
into the two following broad classes |jl|]: 

• Models in which the CKM matrix remains the unique source of both, flavour and 
CP violation. The effects of this source are however modified by the new interaction 
vertices (of new particles) in which the CKM matrix elements appear. 

• Models with entirely new sources of flavour and/or CP violation. 

The first class can be conveniently further subdivided into the so-called MFV models 



2|, I and the generahzed MFV models (GMFV). 

The characteristic feature of the MVF models (MVF scenarios of new physics) is the 
strong dominance in their low energy effective Hamiltonian of the same operators that 
occur in the low energy effective Hamiltonian of the SM. In such models the formula (|1 . 1|) 
remains valid and the relation between the ratio AMg/ AM^ and the length of one side 
of the unitarity triangle, Rt-, is as in the SM: it remains independent of the parameters 
of the particular model. Thus, for this class of models the unitarity triangle is universal 
^. The distinction between different models of this class can then be made through the 
study of e and AM^ which in contrast to a^Ks and AM^/AM^ do depend explicitly on 
new physics contributions. A detailed analysis of the profile of the UT in supersymmetric 
scenarios of this class can be found in [^ . Other discussions of the MFV models can be 
found in 0, gg, |2g . 



The GMFV models generalize the MFV models by allowing for significant contribu- 
tions of the nonstandard operators in the effective low energy Hamiltonian. In this class 
of models the formula ( |1 . 1| ) is still valid but the relation between AM^/AM^; and Rt is 



modified. Hence, Rt determined from tlie measured ratio AMg/AMd does depend on tlie 
parameters of the model. 

While the presence of new CP violating phases in B^ s~^d s mixing and K^-K^ mixing 
could turn out to be necessary to explain the future precise value of a^Ks ; it is important 
to investigate first the scenarios that do not invoke new sources of flavour and/or CP 
violation. It is therefore useful to analyze first MFV and GMVF models that are more 
constrained than the more general scenario mentioned above. 

It turns out that in the MFV models there exists an absolute lower bound on sin 2(3 [|1^ 
that follows from the interplay of AMd and e and depends mainly on | V^b | , | Vub/Vcb \ and the 



non-perturbative parameters Bk-, Fb^\I Bb^ entering the analysis of the unitarity triangle. 
An updated conservative lower bound on sin 2/5 obtained by scanning independently all 
relevant input parameters reads [|ll 

(sin2/3)^in = 0.42. (1.6) 

As analyzed in |1T0|, p|, this bound could be considerably improved when the values of \Vcb\i 
iKife/Kfel, Bki Fb^\I Bb^-i will become better known or if AMg is measured so that the 
ratio AMs/AMrf can be used, along with the non-perturbative parameter S,-, to determine 
the length of one side of the UT. The lower bound ( |1.6| ) is fully consistent with the 
experimental data but as the latter are not yet very precise it could prove useful when the 



knowledge of \Vch\i \VublVch\-, Bk, Fb^x/Bb^, AMg and ^ improves. Note that the bound 



( |1.6| ) allows values of sin 2/? that are slightly smaller than the ones obtained from the fits 
to the unitarity triangle within the SM, (sin2/?)sM > 0.5. In view of the unexpectedly 
high value of a^Xs found by the Belle collaboration [Q, more interesting at present appears 
the upper bound on sin 2/? in MFV models that reads [& 53 



{sm2P)^,, = 2R,^l-Rl. (1.7) 

Here Rb is the lenght of one side of the unitarity triangle (see fig. [^) given in terms 
of iKj^/Kfel in eq. ( p.l2[ ). With the input parameters specified in table 1, one obtains 



(sin2/5)max = 0.82 that is fully consistent with the BaBar result but appears to be 
violated by the Belle result 0. We will return to this issue in the course of this paper. 
The natural next step is to exploit GMFV models. In the present paper we would like 
to make this step and present general formulae relevant for the analysis of the unitarity 
triangle and sin 2(3 in the GMVF models. 

Examples of the MFV models are the Two Higgs Doublet Models II (2HDM(II)) and 
the MSSM, in which sfermion mass squared matrices are aligned with the corresponding 
fermion mass matrices and the CP violating phases of the gaugino masses, fi and A 



parameters are all set to zero, provided the ratio of the vacuum expectation values of 
the two Higgs doublets, ^2/^1 = tan/? is not too large. It is well known |2^, ^ that in 
both these models the contribution of light charged Higgs boson and/or (in the case of 
supersymmetry) charginos and stops to the Wilson coefficient of the standard {V — A) ® 
{V — A) operator can significantly enhance the tW^ contribution to AMg, AM^ and e. 
In this paper we show that for large tan/5 = ^2/^1 both models become GMFV models. 
In particular in the MSSM in the limit M^y < Mjj+ <^ Msparticie, which we consider here 
for simplicity, we find that: 

• There can be significant contributions to AMs from the charged Higgs box diagrams 
(growing like tan^ f3 at the 1-loop level and faster after including leading higher order 
corrections) and, growing as tan^ /3, contributions arising from the double penguin 
diagrams involving the neutral Higgs scalars. Compared to the contribution of the 
extended Higgs and chargino/stop sectors relevant for low tan/?, the interesting 
feature of all these new contributions is their sign which is opposite to the standard 
tW^ box contribution. 

• Compared to AMg, the corresponding contributions to AMd and e are suppressed 
by the quark mass ratios nid/ms and rrid/mb, respectively. 

• Consequently, in this scenario sin 2f3 cannot deviate significantly from its SM value, 
i.e. the lower bound (ll.Gl) can never be reached. 



• Present experimental data strongly limit large mixing of stops if their mass difference 
is large compared to the electroweak scale. 

• If a^Ks is found below 0.5 or above 0.82 this particular supersymmetric scenario will 
be disfavoured (together with the SM). If a^^s ~ 0-''' this scenario can lead to the 7 
angle slightly smaller (depending on the measured value of AMg) than in the SM. 

The full MSSM with large tan/3, including the effects of light sparticles, will be analyzed 
in the forthcoming paper ^ . 



Our paper is organized as follows. In section 2 we generalize the MFV formulae 
of [|l3, ^ m to GMFV models. While in MFV models the new physics short distance 
contributions to i?°g-i?^^ mixing and e can be described by only a single function Fu, 
the transition to the GMFV models (in which new operators contribute) requires the 
introduction of three rea/ functions Ff^, Ff^ and F^f-. We present simple strategies involving 
the ratio AMg/AMd, sin 2/5 and the angle 7 that allow to search for the effects caused by 
new operators (sec. 2.3) and discuss model independent bounds on the functions F^^' '^ 
which follow from the present and future experimental data (sec. 2.4). We also point out 



that in this class of models the function F^^ can be directly measured through AM^ and 
that the angle 7 could be larger than 90°. The formulae necessary to express the functions 
Fif '^ directly in terms of the Wilson coefficients of the four-fermion operators are collected 
in Section 3. The 2HDM(II) and the MSSM with large tan/3 = ^2/^1 and heavy sparticles 
are discussed in Section 4. We give complete formulae for the one loop contribution of 
the box diagrams involving charged Higgs bosons and derive simple approximate formulae 
describing the dominant effects of the double penguin diagrams. Consequences of their 
large contribution and implications of the bounds presented in sec. 2.4 are also discussed. 
We conclude in Section 5. 



2 Basic Formulae 

2.1 Effective Hamiltonian in GMVF Models 

The effective weak Hamiltonian for AF = 2 transitions in the GMVF models can be 
written as follows 

^^P., = ^^ Y. T/^KMa(/i)g. . (2.1) 

Here Qi are AF = 2 operators, Gp is the Fermi constant and Vq-^^-^ is the appropriate 
Cabibbo-Kobayashi-Maskawa (CKM) factor. Because in the GMFV models the CKM 
matrix is by definition the only source of fiavour and CP violation, the Wilson coefficients 
Cj(/i) are real. Using this Hamiltonian with the Wilson coefficients evaluated at the ap- 
propriate scale ^ one can calculate AF = 2 amplitudes, in particular the mass differences 
AMd^s and the CP violation parameter e measured in K -^ titt decays. 

The full set of dimension six operators contributing to AF = 2 transitions consists of 
8 operators. According to the chirality of the quark fields they can be split into 5 separate 
sectors. The operators belonging to the VLL, LR and SLL sectors read: 

Qi"" = {dj-f^PLdi){djYPRdi), 
QLR ^ {djPMidjPRdi), 

gfL = {djPLdi){djPLdi), 

gfL = {dja^,PLdj){dja^''PLdj), (2.2) 

where J, J are the flavour indices (i.e. d^ = b, d2 = s, di = d and, analogously, u^ = t, 
U2 = c, ui = u), a^y = ^[ifMylu], Pl,r = 2(-'- "F 75) ^^^ t^^ colour indices are contracted 



within the brackets. The operators belonging to the VRR and SRR sectors are obtained 
from QY^^ and Qf^^ by interchanging Pl and Pr. Since QCD preserves chirahty, there 
is no mixing between different sectors. Moreover, the QCD evolution factors from high 
energy to low energy scales in the VRR and SRR sectors are the same as in the VLL 
and SLL sectors, respectively. However, one should remember that the initial conditions 
Ci{fit) (where fit = 0{mt)) are in general different for operators involving P^ and Pr. 
In the limit in which the effective Hamiltonian (|2.1| ) is dominated by the single QY^^ 
operator one recovers the results of the MFV models. 

The QCD renormalization group factors relevant for the Hamiltonian ( p.l|) have been 
calculated at the NLO level in ^7\, ^, ^, ^, ^, ^ where the last four papers deal 
with the LR and SLL(SRR) operators. In particular in ref. |^ master formulae for 
AF = 2 NLO QCD factors relating Ci{fit) to Ci{fi) where fit = 0{mt) and yU = 0{mh) or 
fi = 0{2 GeV) have been presented and evaluated numerically in the NDR renormalization 



scheme. Below we will exploit the general formulae of ref. |22[ expressing AM^, AM^ and 
e in terms of the non-perturbative parameters Bi. 
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Figure 1: Unitarity Triangle. 



2.2 AMd, AMs, £ and sin 2/3 in the GMVF Models 



It is straightforward to generalize the formulae of refs. [0, |T^, ||, |2^ to the case of GMFV 
models. To this end, following the notation of ref. 0, we introduce three real functions 



So{xt)[l + fd] 



So{xt)[l + fs 






So{xt)[l + f, 



(2.3) 



relevant for AM^, AMg and e respectively. S'o(xi) with Xt = mf/M^^ is the function 
resulting from box diagrams with (t, W^) exchanges. So{xt) ~ 2.38 ± 0.11 for rht{mt) = 
(166 ± 5) GeV. In order not to complicate the expressions below we assume F^^ > 0. 
Generalization to negative F^^ can be easily done following the discussion of ref. [^ . We 
further split the parameters fi into universal and non-universal parts 



Ji Juni ~r Ji 

where /j = in the MFV models |jlO|- We have then, 



AM„ 



67r^ 



MB^VBBB,F'\Vt,\'Fl q = d,s 



(2.4) 



(2.5) 



where Fsg is the Sg-meson decay constant, Bsg is a non-perturbative parameter and 
riB = 0.55 is the QCD factor [^, |2S[. The measurements of AMq determine the length 
Rt of one side of the unitarity triangle (shown in fig. |1|) defined by 



Rt 



\Vt,Vtl\ 



{l-gy + V 



X 



V, 



td 



Vr 



do 



(2.6) 



Here q = g{l — A^/2), fj = 77(1 — A^/2) ||^, and A, g and r] are the Wolfenstein parameters 
||33|| . As in ref. ^3| we set A = 0.222 in the analytic formulae below.Q Other input 
parameters are collected in table |l|. 
From AMd and AMd/ AMs we find 



Rt = 1.084 



Ru 1 



A 



Rq 



AMd 



0.487/ps 



230 MeV 



^BaFBa 



10.55 
Vb 



and 



Rt = 0.819 ^\ 



AMd 



'15/ps 



0.487/ps V AM, 



R 



sd 



(2.7) 



(2.^ 



*Because of that some numerical factors in the formulae below differ from their counterparts in refs. 
0| where A = 0.220 has been used. This change has only a very small impact on the numerical 



analysis. In particular the bound (1.6) remains unchanged. On the other hand the increased value of A 
shifts Vud closer to its experimental value. See ref. |ll[ for the discussion of this point. 



Table 1: The ranges of the input parameters. 



Quantity 


Central value 


Error 


A 


0.222 


±0.0018 


iKb 


0.041 


±0.002 


Kb/Kfel 


0.085 


±0.018 


IK^I 


0.00349 


±0.00076 


Bk 


0.85 


±0.15 


e 


2.280 X 10-3 


±0.013 X 10-3 


yBB^Fsa 


230 MeV 


±40 MeV 


4^Fb^ 


265 MeV 


±40 MeV 


i 


1.15 


±0.06 


mt 


166 GeV 


±5 GeV 


AAf, 


0.487/ps 


±0.009/ps 


AM, 


> 15.0/ps 




Mw 


80.4 GeV 





respectively where the Wolfenstein parameter A is defined by \Vcb\ = A\^. Here 



Rs 



l + fs 



i 



'Bs 



B 



Bs 



Fb^^/Bb^ 



i + U 

The measurement of the parameter e imposes the constraint which reads: 



(2.9) 



V 



(1 - g)A%Fit + P,{e) A^Bk = 0.204 



(2.10) 



where 772 = 0.57 is the QCD factor [|2^ and Pc{s) summarizes the contributions not 
proportional to Vl*Vtd- With ( |2.7|) and (|2.10| ), the formula of ref. [T^ for sin 2/5 valid in 
MFV models generalizes to 



sin 2/3 



1.65 



-R 



de 



0.204 



A^B 



vPcie) 



K 



R, 



de 



l + fd 
l + fe 



:2.ii) 



Note that new physics can affect sin 2/3 both through fd and fe in (|2.11|) and indirectly 
through fj. We assume as in [|], |10| that new physics contributions to Pci^) are negligible. 
In this case Pc{e) = 0.30 ± 0.05 [0. 



2.3 How to distinguish GMFV from MFV 

In general, the presence of new AF = 2 operators causes the departure from the rela- 
tion Ff^ = F^l = F^^ valid in the MFV models. This means for instance that the ratio 
AMs/AMd, being now dependent on new physics contributions, cannot be used any longer 
for the construction of the universal unitarity triangle ||^. In other words the dictionary 
between Rt and AMs/AMd-, as given by ( |2.8| ), differs from the corresponding one in the 
MFV models because Rsd 7^ 1- This fact offers a possibility to distinguish experimentally 
between these two classes of models. Two strategies are presented below. Because of the 
unitarity of the CKM matrix these strategies are related to each other. 

2.3.1 Strategy A 

For given values of iV^^ib/Kfel and AM^/AMd one can determine (see fig. |l|) 



W V* I A'^ 1 



V. 



ub 



(2.12) 



Vcb 

and Rt by means of ( p.8| ), respectively. This gives the apex of the unitarity triangle with 



and consequently 



e = '^(1 + Rl - Rf) , ri = ^jRl - g'. (2.13) 

sm2/J=M_^. (2.14) 

These formulae establish the relation between sin 2/5 and AMg/AMd that depends on 
the ratio Rsd- In fig. || we show sin 2/3 as a function of AMs/AMd for iKit/Kfel = 
0.070, 0.085, 0.10 and different values of the ratio Rsd- To this end we have set ^ = 1.15. 
We observe that for AMg/AMd below 40 the distinction between GMFV and MFV models 
will be difficult unless the values of C,, \Vub/Vcb\, a^Ks and AMg will be known very accu- 
rately or the value of Rsd differs substantially from unity. For larger values of AMs/AMd 
the distinction is clearer. 

From fig. ^ the impact of the measurement of AMg/AMd on the allowed values of 
sin 2/5 is clearly seen. For 0.8 < Rgd < 1-2 and AMs/AMd < 40 only values compatible 
with current experimental result ( |1.3| ) and well above the bound ( |1.6| ) are allowed. On 
the other hand, for sufficiently low or sufficiently large values of Rsd, smaller values of 
sin 2/? are also possible. While such low values of sin 2/? are still compatible with the 
BaBar result, they seem to be excluded by the Belle measurement of a^Ks- In fact, as 
follows from fig. H, the latter measurement can be accomodated only if Rsd ~ 1 and 
\Vub/Vcb\ ~ 0.1. 
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Figure 2: sin2/5 as a function of AM./AM^ for |Kfe/Kb| = 0.070 (panels a and d), 0.085 
(panels b and e), and 0.10 (panels c and f) for different values of Rgd (marked on the 
curves) and C, = 1.15. 

It is also easy to find that the two possibilities, small Rsd and large Rsd, favour small 
and large values of the angle 7 in the unitarity triangle, respectively. Which of these two 
possibilities is favoured by the data can only be decided by other measurements. This 
includes e, AM^ alone and in particular a direct measurement of 7. This brings us to the 
second strategy. 
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2.3.2 Strategy B 

The angle 7 in the unitarity triangle can be found from Rt and /3 determined through 

( p. . 1| ) by using the relation 

1- RtCos/3 

cot 7= p ■ ^ • 2.15 

Rt sm p 



Expressing Rt in terms of AMs/AMd by means of ( |2.8|) allows to calculate 7 as a function 
of sin 2/3, AMs/AMd and Rsd- For i^^^ ^ 1 the predictions for 7 in GMFV models will 
generally differ from those in the MFV models. Comparing these predictions with future 
direct measurements of 7 it will be possible to distinguish between these two classes of 
models and check whether the inclusion of new operators is required by the data. In fig. |] 
we show 7 as a function of AMg/AMd for sin 2/3 = 0.4, 0.6, 0.8 and various values of 

Rsd- 

We observe that the distinction between MFV and GMFV models in this strategy is, in 
contrast to strategy A, very transparent in the full range of AM^/AM^ considered. As this 
strategy involves only a^Xs and AMs/AMd, that are theoretically cleaner than |\4fe/Vcb|, 
it is this strategy which in the future should play the crucial role in the distinction between 
the MFV and GMFV models. The ratio AMs/AMd and the asymmetry a^xs should be 
determined very precisely in the coming years. The determination of 7 is more difficult 
but should be achieved at LHCb and BTeV. Some information on the angle 7 should also 
be gained from the Bd — > ttK decays measured by CLEO, BaBar and Belle and by the 
combination of Bd -^ tt+tt^ rate (already measured by these three collaborations) and 
the rate of the Bs — > K^K' decay [^ which are going to be measured at Tevatron. 

As seen in fig. ^ the values of 7 for Rsd < 1.2 are below 90°. On the other hand for 
substantially higher Rsd also 7 > 90° is possible. The possibility of 7 > 90° resulting from 



the unitarity triangle fits is very interesting in view of several analyses ||36| , p7| , pSJ p9|] 
of two-body non-leptonic decays B —>■ ttK, tctt that favour 7 > 90° in contradiction 
with the usual unitarity triangle analyses that confidently give 7 < 90°. With increasing 
AMs/AMd this problem will become more serious. 

In view of sizable theoretical uncertainties in the analyses oi B ^ ttK, vrvr and of large 
experimental errors in the corresponding branching ratios it is not yet clear whether the 
discrepancy in question is serious. For instance ||4^ sizable contributions of the so-called 



charming penguins to the B — > nK amplitudes could shift 7 extracted from these decays 
below 90° but at present these contributions cannot be calculated reliably. Similar role 
could be played by annihilation contributions and large non-factorizable SU{3) breaking 
effects ||3B|. Also, a new physics contribution in the electroweak penguin sector could also 
shift 7 to the first quadrant [0 . It should be however emphasized that the problem with 
the angle 7, if it persisted, would put into difficulties not only the SM but also the full 
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Figure 3: 7 as a function of AMs/AM^ for sin 2/? = 0.4, 0.6, 0.8 and various values of 
Rsd (marked on the curves). 

class of MFV models in which the lower bound on AMg/AMd implies 7 < 90°. On the 
other hand as seen in fig. ^ for sufficiently high values of Rsd, the angle 7 resulting from 
the unitarity triangle analysis can easily be in the second quadrant provided AMg/AMd 
is not too large. 

Clearly a general analysis of the unitarity triangle involving e, AMg^d, |Kife/Kb| and 
|Vcfe| can also be used to search for the effects caused by the new operators but the two 
strategies outlined above have in our opinion the best chance to distinguish between MFV 
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and GMFV models in a transparent manner. 



2.3.3 Unitarity Triangle, sin 2/3 and 7 



A different version of strategy B is to construct the unitarity triangle by means of the 
ratio AMs/AMd and the asymmetry a^Ks- With Rt given by ( p.8|) the parameters g and 
f] can be determined from the formulae[| 



1 + sin 2p + y^l - sin 2p] (2.16) 



= 1 

'^ 2 

(recall that in in GMFV models sin 2/? = a^Ks) obtained directly from eqs. ( |2.6| ) and 
( p.l4| ). These formulae are equivalent to those presented in [|], ^ but are more elegant. 

As an illustration we show in fig. ^ the ranges of {g, ff) allowed by the hypothetical 
measurement AM^ = (18.0 ± 0.5)/ps for three values of a^^s ^^^d different values of Rsd- 
Solid ellipses correspond to Rsd = 1 valid in particular in MFV models. We also show the 
i?fc-constraint, eq. ( |2.12| ), with i?{, = 0.37 ± 0.08 and as a useful reference the e-constraint 
( p.lOp with 1 + /e = 1 corresponding to the SM. 

From fig. ^ it is clear that for a^Ks = 0.80±0.05 the class of models giving Rsd = 1 and 
1 + fe = 1 (which includes also the MFV models) is consistent with all constraints but only 
for sin 2f3 in the lower part of the chosen range. There is also a room for contributions 
of new operators resulting in Rsd 7^ 1 provided 0.7 < Rsd < 1-4. As in strategy B 
they could be distinguished through the value of the angle 7. We also observe that no 
sizable contributions to e beyond the SM ones are required. For a^^s = 0-60 ± 0.05 
models giving Rsd = 1 are consistent with the Rb- and e-constraints for sin 2/5 in the 
full chosen range. In this scenario models with 0.8 < Rsd < 2.0 and no sizable new 
contributions to e and models with Rsd < 0.8 but with 1 + /^ > 1.0 are favoured. Again 
the measurement of 7 could distinguish between these possibilities. It is interesting to 
note that for 1.5 < Rsd < 2.0 it is possible to have 7 > 90 even for 1 + /^ = 1. Finally, for 
atpKs = 0.4 models giving Rsd = 1 are ruled out as they do not satisfy the i?b-constraint. 
In order to reconcile this constraint with a^Ks ~ 0-4, values of Rsd substantially different 
from unity are required. Moreover in the case of Rsd < 1-0 one has 7 <^ 90° but large 
new contributions to e leading to 1 + /^ > 1.0 and Rde < 1-0 are mandatory. In contrast, 
if Rsd ^ 2.0 7 can be much bigger than 90° even without new contributions to e i.e. with 
1 + /. ~ 1. 

tThis is valid for — 7r/4 < (3 < 7r/4; for other ranges of /? similar formulae can be obtained. 
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Figure 4: Ranges of {p,f]) allowed in la for AMg = (18.0 ± 0.5)/ps, three values of a^Xs 
and different values of Rgd (marked in the figures). Black spots correspond to Rgd = 1- 
Dotted lines show the constraint from the parameter e, eq. ( p^.lO| ), for 1 + f^ = 1. 
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Figure 5: Lower bound on sin 2/3 in GMFV models as a function of the ratio Rde- Dashed 
hues show the l-sigma BaBar result (sin 2/5 = 0.59 ±0.15), the horizontal solid line shows 
the lower limit from the Belle measurement (sin 2/5 = 0.99 ± 0.15) and the shaded area 
corresponds to the sin 2/? range presently allowed in GMFV models at la by the official 
grand average ( |1.3D . 



As in some scenarios discussed above Rde must differ from unity, we show in fig. ^ the 
dependence of the lower bound for sin 2/5 on the value of Rde together with the la BaBar 
result and the present official experimental la band ( |1.3| ) for a^Xs = sin 2/?. It follows, 
that at present the ratio Rde would be constrained at la by the BaBar result to be less 
than 2.2 but the analogous limit following from the grand average (|1.3|) is much higher. 
Fig. 1^ has been obtained for Rgj, = 1.0 but (sin2/?)inin depends only very weakly on the 
ratio Rsd'- very similar curves are obtained also for 0.6 < Rgd < 2.0. 



2.4 Constraints on GMFV models from unitarity triangle 

In the preceding subsection we have presented two strategies which in principle should 
allow to decide on the basis of experimental measurements whether going beyond the 
MFV models is necessary, i.e. to establish whether Rgd 7^ 1- In this section we want 
to explore constraints and correlations imposed by the experimental data (present and 
future) and the unitarity of the CKM matrix on the functions F^^. These constraints can 
be then effectively used to test the specific GMFV models of new physics. 

The first constraint follows from fitting the formula ( p.5|) to the measured (in the near 
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future) value of AMg. This determines 1 + /« 



(or F/J 



1 + /, = 0.80 



2.38 

So{xt) 



265 MeV 



BbsFbs 



2 


"0.55' 




"0.041' 

. Vts\ _ 


2 


"AM/ 
_15/ps_ 



(2.17) 



This formula follows also by equating Rt determined from AM^ alone (eq. (|2.7] )) and 
Rt determined from the ratio AM^/AM^ (eq. (|2.8| )). Scanning over uncertainties in the 
ranges specified in table ^ and setting \Vts\ = \Vcb\ gives 



0.52 



AM/ 
15/ps 



<l + fs< 1.29 



AM, 
15/ps 



(2.18) 



(At present this gives of course 1 + /, > 0.52.) It is worth emphasizing that this bound 
is independent of the uncertainties of iV^b/Kbl as well as of any assumptions about pos- 
sible new physics contribution in the K^-K^ system. Therefore the GMFV models not 
respecting it are (will be) ruled out. We will see in sec. 4.2 that the MSSM, for some 
values of its parameters, violates precisely this bound. 

Next, there are bounds on Rt coming from the requirement that 



l-Rb<Rt<l + Rb 



:2.i9) 



which gives 0.54 < Rt < 1.46. This can be used to constrain either 1 + fd or Rsd depending 
on how one determines Rt. In the first case one gets 



0.20 <l + fd< 4.24 



and in the second 



0.29 



AM, 
15/ps 



< Rsd < 2.73 



AM, 
15/ps 



(2.20) 



(2.21) 



More stringent constraint on Rsd (if Rt is determined from ( |2.8D ) or 1 + /^ (if Rt is 
determined from (|2.7| )) will follow from Ri, combined with the information about sin 2/5 
obtained from future accurate measurements of the asymmetry a^^s- I^ is easy to see 
that for sin 2/5 < 0.34 there are two allowed bands of 1 + /^ corresponding to two possible 
solutionsfl for Rt. 

R^ = cos /3 T ^Rl - sin^ p (2.22) 

(the solutions with — (+) correspond to smaller (larger) values of the angle 7). For larger 
values of sin 2/? the two bands overlap which means that there is only one allowed range of 



mm 

b ■ 



■l-The band of allowed Rt values splits into two for sin 2/? « 0.56 corresponding to sin/3 = R 
Additional uncertainties in translating Rt into 1 + fd (or Rsd) result in lowering the value of sin 2/3 below 
which the two allowed bands of 1 + /d (or Rsd) appear. 
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Figure 6: Allowed bands of Rsd as a function of AM^/AM^ for different values of sin 2/?. 



1 + /d (or Rsd)- For example scanning over uncertainties one obtains: 0.21 < 1 + fa < 0.78 
or 0.84 < 1 + fd < 4.15 for sin 2/? = 0.2, 0.23 < I + fd < 3.85 for sin 2/? = 0.4, 
0.26 < 1 + fd < 3.27 for sin 2/3 = 0.6 and 0.44 < I + fd < 1.98 for sin 2/? = 0.8. The 
corresponding allowed ranges of Rsd are shown in fig. |^ for different values of sin 2j3 as 
functions of the measured values of AMs/AMd. 

Further constraints correlating 1 + fe with 1 + fd can be obtained by using the exper- 
imental information about the parameter e. To this end, with a range of Rt determined 
from 1 + fd and AMd by scanning over the relevant uncertainties, one checks whether 
there exist values of p and fj satisfying eqs. ( p.6|) and (|2.10|) for a given value of 1 + /e. It 
is easy to see that from AMd and e alone only a very weak lower bound on |1 + /e| can 
exist for 1 + /^ ~ 1.8. As a next step, one can impose the constraint from Rt, (eq. ( |2.12| )). 
The resulting allowed range in the plane (l + Z^^, 1 + fe) is shown in fig. ^,b by the dashed 
lines (their vertical parts correspond to the lower bound ( |2.20| )). For 1 + /^ in the range 
( p.20| ) no upper limit on |1 + /^l from AMd, £ and Rb exists (except for a very narrow 
range 0.7 < 1 + fd < 0.85). This is because for almost all values oi 1 + fd satisfying (|2.20|) 
the range of possible values of Rt is such that it is possible to satisfy the constraint from 
Rb with f] = which in turn allows to suppress arbitrarily large values of 1 1 + /^ | and to 
satisfy the eq. ( p.lO|) . 

It is also interesting to assume that the asymmetry a^^s (i-^- sin 2f3 in GMFV models) 
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Figure 7: Allowed ranges of 1 + /^ and l + fe- AM^, e and Rb allow the region delimited by 
the dashed lines. Regions between the solid lines are allowed by AM^, e and sin 2/5 = 0.4 
(panel a) and sin 2/5 = 0.8 (panel b). Dotted regions are allowed by AM^, e and Rj, for 
sin 2/? = 0.4 and 0.8 in panels a) and b), respectively. 



is measured with sufficient accuracy and to correlate l + fe with 1 + /^ by using the 
experimental information about AM^ and e for fixed values of sin 2(3. To this end, for Rt 
given by eq. (|2.7| ) one determines the parameters p and f/ from the formulae (|2.16|) and 
checks whether eq. (|2.10| ) can be satisfied. Constraints on the plane (1+/^, l+/e) obtained 
in this way are shown in fig. ^a for sin 2/5 = 0.4 and in fig. |^ for sin 2/5 = 0.8 by the 
solid lines. It should be stressed that it is the constraint from the a^xs asymmetry which 
eliminates solutions with f] ^ thus providing, for fixed 1 + fd, the upper bound on l + Z^. 



In agreement with the bound ([O]) we observe that for sin 2/5 = 0.4, Rde = {l + fd)/{l + fe) 
has to deviate from unity while for sin 2/5 = 0.8 points corresponding to Rde = 1 lie within 
the allowed region. 

Finally, one can impose also the constraints from Rb ( p.l2| ). The allowed ranges in 
the (1 + fd, 1 + fe) planes are dotted in figs. ^,b. They are not simply given by the 
intersection of the regions allowed respectively by (AM^, e, Rb) and (AM^, e, sin 2/5) 
because the same point in the (1 + fd, 1 + fe) plane may require different p and f/ to be 
compatible with the two above sets of experimental data. 

For a fixed value of the ratio AMg/AMd the same analysis can be also repeated for the 
{Rsd, 1 + fe) plane. However, since the value of AMg/AMd serves only to determine Rt 
from eq. (|2.8|) we show instead in fig. |^ the allowed ranges in the plane {nRg^, l + fe) where 
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Figure 8: Allowed ranges of nRsd and 1 + /e where /t = xirfAAr- '^^s/^Ma, e and 
i?f, allow the region delimited by the dashed lines. Regions between the solid lines are 
allowed by AMs/AMa, e and sin 2/3 = 0.4 (panel a), sin 2/3 = 0.6 (panel b), sin 2/3 = 0.8 
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K 



30 



aMs/am'- ^^ ^^^^ figure, panels a)-c) correspond to sin 2/3 = 0.4, 0.6, 0.8 respectively. 
Panel d) corresponds to the averaged experimental value 0.80 ± 0.11. Because of the 
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absolute constraint ( ^.IT] ), Rsd oc 1/(1 + /d) and qualitatively figs. ^ and |b are obtained 
from figs. 0a and ^, respectively by subjecting the x-axis to the transformation x -^ 1/x. 
Quantitatively however, the bounds on GMFV models provided by fig. ^ are tighter as 
they are independent of the poorly known parameter Fb^vBbj^- 



General formulae for F^, 
Models 



F^ 



and F/^ in GMFV 



Using general formulae (7.27)-(7.32) in [^] it is easy to express the functions Ff^, F/^ and 
F^^ in terms of the Wilson coefficients at the scale at which the effective Hamiltonian is 
generated, the relevant QCD renormalization group factors r] and the non-perturbative 
Bi factors. Suppressing the superscripts d, s and e for a moment we find 

1 



VLL 



So{Xt) + T-C'new 



(/^t 



4^; 

jSLL 



VRR/ 



LR/^LR/ 



(/.,) + prc^r(/^*) + ^2"^c72 ^(/^ 



SRR/ 



^p«LL Cf^^^^,)+Cf^^{^^,) 



Pi 



SLL 



^SLL 



(^t 



+ C^f ^(/i.^ 



(3.1) 



where r = 0.985 |2^ describes the QCD corrections to So{xt) in the SM. This factor is 
present because we have factored out rjB in (|2.5| ) and ri2 in the analogous formula for e. 
The first line of eq. ( |3.1|) contributes to funi in eq. (|2^ ) and is therefore present also 
in the MFV models. The remaining lines are characteristic for the GMFV models and 
contribute to the /, in eq. ( p.4|) . Thus, different values of the functions Ff^, F^^ and F^^ 
originate from generally different values of the Wilson coefficients C°-{fit) and factors P" 
pertinent to AM^, AMs and e: 



pa 



P^/{Ai^bBb, 

Pt/{^V2BK) 



{AMd,s 



(3.2) 



where Bb^ ^ and Bk are the relevant non-perturbative parameters related to the matrix 
elements of Q^^^- 

In the case of F^i'^ the coefficients P^ are given by 



pLR 



P 



LR 



^l' 



SLL 



'\ [^ll(/^fe)]LR [B¥'{^^b)\^^ + \ hl(/W6)]LR [^2''(/^fc)],ff ' 


(3.3) 


"2 [^12(/"fc)]LR [BW^^l>)\^^ + 4 [^22(/i6)]LR [^2 ""(/^b)] ,ff > 


(3.4) 


^ [^ll(/^fe)]sLL [^l'^'^(/^t)]eff - 2 [^2l(/^6)]sLL [^'"^'^(Z^'')] eff ' 


(3.5) 
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p. 



SLL 



I [^12(/^b)]sLL f^f '^(/"fe)l„„ - I [^22(/i6)]sLL \Bf''{^^b 



off 2 
with the effective parameters [-Bf(/^b)]efr defined by 



eff 



(3.6) 



[Bt{^^b)U = 



M, 



Ba 



1.44 



\mb{iJ,b) +mq{iJ,b) 

4.4 GeV 
nibifib) +mg{^b) 



Btifib 



2 r 



M 



n2 



B« 



5.28 GeV 



i?f(;"6 



(3.7) 



where Bf{iib) are related to the hadronic matrix elements {B^\Qi\B^). The QCD factors 
[Vij{^^b)]a are given in [||]. 

In the case of Ff^ the QCD factors ["riij{l^b)]a ^^^ replaced with [rjij{^L)\a where /i^ = 
2 GeV and the corresponding effective parameters [B^^iil)]^^ are defined by 



[B■{^^L)] 



rriK 



cff 



msifXL) +md{iJiL] 



Bti^L 



18.75 



115 MeV 



n2 



rUsifiL) +md{lJ,i 



BK^l). (3.8) 



The NLO QCD factors fjijifJ^b) and 77jj(2 GeV), relevant for AM^^s and e respectively, 
have been evaluated in [^ . For completeness we give in table ^ their numerical values 
for different scales (denoted here by /xnp) at which the new operators are generated. 

The parameters Bi for the K^-K^ mixing are known from lattice calculations [^l], |41 
For the values corresponding to ^t in the table El and in the NDR scheme one finds ||32 



LR 



Pi 

pSLL 



-36.1, 
-18.1, 



pLR 
-'2 

pSLL _ 
-'2 




for 



/i = 2 GeV. 



(3.9) 



The large values of these coefficients originate in the strong enhancement of the QCD 
factors rjij for the LR and SLL (SRR) operators and in the chiral enhancement of their 
matrix elements seen in eq. ( p.8|) . Consequently even small new physics contributions to 
Cj"^{fit) and Cf^^{fit) can play an important role in the phenomenology pT| , |i2 . 

In the case of the B^-B^ mixing the chiral enhancement of the hadronic matrix ele- 
ments of the LR and SLL operators is absent. Moreover, the QCD factors rjij are smaller 
than in the case of the K^-K^ mixing. Consequently the coefficients P}"^ and pf^^ are 
smaller in this case but can be still important. As lattice results are not yet available for 
the hadronic matrix elements of the LR and SLL operators in the B system |^3| we will 
set in this case B^i^Hb) = 1. Taking Mb = 5.28 GeV, fib = 4.2 GeV, mb{fib) +fnq{fib) = 4.2 
GeV and as{Mz) = 0.118 one finds 



pLR ^ _x 55^ 
pSLL ^ _^_49^ 



pSLL __ 




for 



fJ-b 



4.2 GeV. 



(3.10) 



21 





/ife = 4.2 GeV 


IJL = 2 GeV 


AiNP 


f^t 


500 GeV 


1 TeV 


f^t 


500 GeV 


1 TeV 


[^]VLL 


0.838 


0.809 


0.793 


0.787 


0.759 


0.744 


[^ii]lr 


0.919 


0.907 


0.902 


0.906 


0.900 


0.898 


[^12] LR 


-0.043 


-0.054 


-0.060 


-0.089 


-0.107 


-0.118 


[^2i]lr 


-0.919 


-1.190 


-1.360 


-1.548 


-1.923 


-2.159 


[^22] LR 


2.303 


2.701 


2.951 


3.227 


3.785 


4.136 


[^lllsLL 


1.676 


1.846 


1.949 


2.063 


2.272 


2.398 


[^12]sLL 


2.049 


2.470 


2.715 


2.970 


3.441 


3.717 


[^2l]sLL 


-0.007 


-0.008 


-0.009 


-0.009 


-0.011 


-0.012 


[^22]sLL 


0.540 


0.480 


0.449 


0.414 


0.366 


0.341 



Table 2: Numerical values for the ?7-factors for the B -B and K -K mixing for 
af\Mz) = 0.118 and different values of the scale /xnp at which the New Physics is 
integrated out. 



Since we took fih = 4.2 GeV instead of 4.4 GeV, these numbers differ slightly from those 



given in ref. ||3^. Finally in order to calculate Pf in ( |3.2| ) we will use rjE = 0.55, ■r]2 = 0.57 

Ex = 0.85 ±0.15, 5b, ^ = 1.30 ±0.18 . (3.11) 
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4 Ff^, F^^ and F/^ in realistic GMVF Models 

4.1 2HDM(II) with large tan/5 

To see what values of Ff^, F^^ and F^^ can be realized in realistic GMFV models we consider 
here two extensions of the SM. The first is the 2HDM(II) which introduces three neutral 
physical scalars (/i°, H^ and A^) and a charged physical scalars H^. At one loop only 
the charged scalars are relevant for the box diagrams contributing to K^-K^ and B^-B^ 
mixing amplitudes. Using the compact notation of ref. |^4| the tree level couplings of 



H^ = {H~^, G^) (where G^ is the would-be Goldstone boson) read 



Alk ; 



Lint = H+UAVAi{ar'PL + a^"'PR)di + He. 



(4.1) 
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where 



-,AIk 



e m 



UA 



V2swMw 



X 



cot /5 

1 



for k = 1 
for k = 2 



(4.2) 



Alk 



e rudj 



V2swM^ 



X 



w 



tan/? for k = 1 
- 1 for A; = 2 



(4.3) 



(Recall that in our notation ds = b, d2 = s, di = d and, analogously, u^ = t, U2 = c, 
Ui = u.) Contribution of H^ to the Wilson coefficients Ci of the operators responsible for 
the transition djdj -^ djdj in eq. ( f2.1| ) can be easily expressed in terms of the coefficients 
a^^^ and a"^^^. Diagrams with one W^ and one H^ giveQ 



tJk tlh 



G|M^5(+)C2^^(/XNP) 



0„2 
^2 2 



E 



'vy fc=i 



D2{Mw,Mjj+,mt,mt) 



-2D2{Mw, Mjj+,mt, 0) + D2(Mh/, M^+, 0, 0) (4.4) 

k k J 

where the four-point functions Dq and D2 are defined in the Appendix. Diagrams with 
two H^ givep]: 

G|M^5(+)Ci^LL(/iNp) = ^E«i''aL'a*/'aL'^2(M^+,M^+,m,,m,) 



A:,/ 
2 



G|M^5Wcr^(/iNP) = ^E 



tjltlktjktil 
"ij "ij "i? "_R 



fc,Z 



D2{Mfj+,Mfj+,mt, rrtt) 



-2Z}2(M^+, M^+, m^, 0) + D2(M^+, M^+, 0, 0) 
G'|M^5Wc{^^(/iNp) = ^E«f «LV/*^4'A(M^+,M^+,mi,m,) 

G'|M^<5(+)CfLL(/iNp) = ^E«R'4''«*#4V^o(M^+,M^+,mi,mi) 

G^M^^Wcf ^(/iNp) = ^E«*/'4'«*/'4W^o(M^+,M^+,mi,mO 

G2,m2,<5(+)CLI^(/xnp) = E «*/'«L'a*/'=aM'm2Do(M^+, M^+, mi, m,) 

' ^ I k 

kl 



(4.5) 



§The contribution of G* to C}^'^'^ is already taken into accomit in the function 6*0 (a;t). Masses of the 
u and c quarks are neglected. 

^In the sum over k and I in the expression for S'^^^'C^^^ the contribution of G^G^ is excluded. It is 
taken into account in the function So{xt). 
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At one loop there are no contributions to the Wilson coefficients of the tensor operators 
Q^^ and (52^^- In the computations we take /inp = Mh+ and apply the formulae 



given in Appendix C of ref. ||32|. As our calculation of the Wilson coefficients at /xnp 
does not include 0{as) corrections and the relevant matrix elements of LR and SLL 
operators in the 5-system have been evaluated using the vacuum insertion method, there 
are inevitably unphysical scale and renormalization scheme dependences present in our 
final results. We expect that these dependences are small at scales C(/iNp) as the strong 
coupling as(/iNp) is small. They could turn out to be more important at /i = /i;, where a^ 
is bigger. Consequently the evaluation of the hadronic matrix elements of the LR and SLL 
operators relevant for B^-B^ mixing in the NDR scheme is very desirable. This would not 
only remove the unphysical dependences in question but would also give the actual values 
of the relevant matrix elements in QCD. Still we believe that our calculation captures the 
correct size of the dominant new physics effects. 

For large tan/3 and Mh+ k. rrit the leading terms of the above contributions to the 
Wilson coefficients Q are of the order {e'^/{32s^M^) = GlM^): 

3 3 mf 



H 



r\j 



for diagrams with W H^, and 

5Wcf^L~0, ^Wcf^-O, 5WC2^^~-^^^^^tan2/3 (4.7) 

for diagrams with H^H^ . It is clear that for large tan/? the biggest contribution appears 
in (5*^^^C^^. It is of the opposite sign than the contribution of the tW^ box diagram and 
can be significant only for the B^-B^^ transition amplitude for which it is of the order 

where we have used mh{^t) ^ 3 GeV and ms{^t) ~ 61 MeV. Compared to the esti- 
mate of eq. (|4.8|) , similar contributions to 5^^^C2^ for B'^-B^ and K^-K^ transitions are 
suppressed by factors md/ms and md/mb, respectively. As discussed in Sec. 3 all these 
contributions are further enhanced compared to the standard ones by the QCD renormal- 
ization effects and, in the case of the K^-K^ transition, also by the chiral enhancement 
of the corresponding matrix element. As a result, in the case of B^-B^ mixing the con- 
tribution of the (^2^ operator can compete with the contribution of the standard Q^^^ 
one for light charged Higgs boson and large values of tan /?. To demonstrate it we plot in 
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fig. 1^ the value of 1 + /s in 2HDM(II) for different values of its parameters. Fig. ^a shows 
that a significant decrease of 1 + /^ below unity is possible only for Mh+ ~ 100 — 200 
GeV and tan/3 close to its upper limit following from the requirement of perturbativity 
of the bottom quark Yukawa coupling. The corresponding effects in 1 + /^ and 1 + /e are 
negligible. The increase of 1 + /^ above unity seen in fig. ^ for tan/5 < 10 reflects a well 
known universal contribution of the box diagrams to the Wilson coefficient of the QX^^ 
operator which gives l + /s~l + /d~l + /e>l. 
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Figure 9: 1 + /s in the 2HDM(II): a) as a function of tan/5 for Mh+ = (from below) 150, 
250, 300 and 350 GeV and b) as a function of Mh+ for tan/5 = (from above) 40, 60, 80 
and 100. 



Unfortunately, recent refinements in the computation of the b ^ s'y rate [^] together 
with the new CLEO experimental result for this process |^ BR{B -^ ^s7) = (3.03 ± 
0.40 ± 0.26) X 10"^ set the bound Mh+ > 380 GeV g^, ||. This means that in the 
2HDM(II) for the still allowed range of charged Higgs boson masses the decrease of 1 + /^ 
can be very small. Consequently, the SM analysis of the unitarity triangle based on e, 
AMd and AMg is practically unchanged in the 2HDM(II) for large tan/5 ^ 50. 

However, the bound on Mh+ from the b —>■ sj rate does not apply in the MSSM which 
we consider in the next subsection. 
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4.2 MSSM with large tan/3 

As a second realistic GMFV model we consider the MSSM. At the one loop level the 
contributions to the Wilson coefficients of the |AF| = 2 operators ( p.2| ) in the MSSM 
are given[] by chargino-top squark box diagrams and by box diagrams with the charged 
Higgs boson. Since the Higgs sector of the MSSM is (at the tree level) a special case of 
the general 2HDM(II) considered in subsection 4.1, the latter contribution is described 



by the formulae ([4.4|) , ( [4.5|) . It is well known |2J, g^, ^ that for tan/3 not too big, the 
MSSM is of the MFV type and both, the chargino-stop and the Higgs sectors give positive 
contributions to 1 + /^ ~ 1 + /^ ~ 1 + /^ which are the bigger the lighter are particles 
in the loops and the smaller is the value of tan/3. In this section we want to consider 
the MSSM with large values of tan/? which are favoured both by the LEP limit on the 
mass of the lighter neutral Higgs boson and by the recent measurement of the anomalous 
magnetic moment of the muon ^^. As the treatment of the full MSSM contribution is 
complicated, we will for simplicity consider here only the limit of heavy sparticles and will 
concentrate only on the most spectacular effects. Complete analysis of the MSSM will be 
presented elsewhere ^6 . 



In the hmit of heavy sparticles (which is practically realized already for Msparticies ~ 500 
GeV) the one loop diagrams involving charginos and stops are negligible. It is however 
known that for large tan /3 even if sparticles are heavy they can still compensate the H^ 
contribution to the b ^ S'j amplitude allowing for the existence of a light, ~ (9(150 GeV), 
charged Higgs boson |^, |51|. From fig. ^ it follows therefore that, even forQ tan/5 < 50 
and already at the one loop level the contribution of the MSSM Higgs sector to the C^^ 
Wilson coefficient can be non-negligible. 

At the two loop level one has to take into account not only the 0{as) corrections 
to the Wilson coefficients (which in the MSSM arise from exchanges of gluons as well as 
gluinos) but also the dominant two loop electroweak corrections (proportional to large top 
and, in the case of large tan/5, bottom Yukawa couplings). The gluonic 0{as) corrections 
to the charged Higgs box diagrams are expected to be of the same order of magnitude 
as the gluonic correction to the SM t — W^ box diagrams, i.e. moderate |^. Also most 



of the two loop diagrams involving sparticles will give contributions suppressed by the 
inverse of the large sparticle masses. The one loop effects of the MSSM Higgs sector are 
however enhanced by an important class of two loop corrections involving sparticles. In 
the limit of heavy sparticles these corrections can be most easily identified in the effective 



II In accordance with the general framework of this paper we assume here that the CKM matrix is the 
only source of flavour and CP violation. 

**In the MSSM one usually constrains tan/3 to be less than 50 — 55 by requiring perturbativity of the 
Yukawa couplings up to the GUT scale ^ 10^^ GeV. 
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Lagrangian approach [^. While direct contribution of heavy sparticles to the Wilson 
coefficients can be neglected, it is well known that for large tan/^ such sparticles do not 
decouple entirely [Q. Integrating them out modifies the original couplings of the two 
Higgs doublets to the known fermions leading to the following quark Yukawa interactions: 

- e,,Hfu\Y^'qj, - h\'^* u\{^,Y^)^' qn + H.c. (4.9) 

where H\ and H^ are the two Higgs doublets giving at the tree level masses to the 
down- and up- type quarks, respectively and gj, df' and u'^ are the fields of the left handed 
fermions (for simplicity we use here the Weyl spinors). Dominant corrections {^yXdY"^ 
and ( ArfYu)^^ are finite and calculable in terms of the sparticle parameters. In the effective 



Lagrangian approach they can be obtained from diagrams shown in fig. [1^ in which Q, 
V^ and D'^ are the scalar superpartners of the SM fermions g^, d'^ and u'^, g is the gluino 
and hI and H^ are the fermionic superpartners of the two Higgs doublets. 
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Figure 10: Diagrams giving rise to dominant, tan/3 enhanced, corrections A^lrf (diagrams 
a and b) and A^y^ (diagrams c and d). 

The effects of the corrections Ay are threefold. Firstly, they modify the 2HDM(II) 
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relations between the masses mdji^Ht) and the eigenvalues of the Yukawa matrices Y^ 



I A. 






where 6^(1) can be found in ref. |S^. Secondly, they induce additional, oc tan/3, terms 



in the couplings (^^), ( |4.3| ) [p^ , |50| . Both these corrections should be taken into account 
in vertices of the charged Higgs boson box diagrams and constitute, therefore, the two 
loop corrections to the Wilson coefficients of the |AF| = 2 operators. It turns out, 
however, that the most important (for non-negligible mixing of the top squarks) effect of 
the corrections AF is the generation of the flavour non-diagonal, tan/3 enhanced couplings 



of the neutral Higgs bosons to down-type quarks p4| , p5| . In the effective Lagrangian 
approach these couplings originate from the diagram ITHb. Details of their calculation have 



been presented in ref. [^]. They can be also computed diagramatically as in [^6|, ^ what 
allows to take fully into account the complicated composition of charginos. Additional 
contributions to the Wilson coefficients Cf^^, Cf^^ and C^^ are then generated by the 



double penguin diagrams shown in fig. |TT| in which the neutral Higgs bosons are exchanged 
at the tree level between two effective flavour changing vertices generated at one loop. 

The single neutral Higgs penguin diagram for the ((i/)L(R) -^ ('^j)r(l) transition (where 
L and R refer to the quark chiralities) grows as tan^/5 and is proportional to nidj {fndj) 



m:, 



| 56| . Consequently the double penguin diagrams in flg. |TT| grow like tan'^/^ times //t^^, 
rn^ and mdjmdj, respectively. This mismatch of powers in tan/3 and powers of light 
quark masses in Higgs penguin diagrams should be contrasted with box diagrams, where 
each tan/? is accompanied by a light quark mass as seen in ( |4.(j| )-( ^^ . In this manner 



the two-loop electroweak double Higgs penguins can potentially compete with the one 
loop electroweak box diagrams. In the standard diagramatic approach to the calculation 
of the flavour changing djdiH^{h^,A^) vertex the mismatch of powers in question can 
be understood simplest as follows. The diagrams for the genuine 1-PI vertex corrections 
(with physical squarks and charginos in the loops) contribute only terms proportional to 
one power of tan/3. They are however accompanied by two tree level flavour conserving 
djdj-B.iggs and didj-Higgs vertices with flavour violating self-energies (involving squarks 
and charginos) on the external quarks lines. As the external momenta can be neglected, 
the internal fermion propagators of dj and dj cancel respectively the rridj and mdj factors 
present in the quark-quark-Higgs vertices. 

For the transitions djdj —>■ djdj the dominant terms obtained from the double penguin 
contributions are 
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(4.11) 



OLEM m; 



2nsL M, 



-md,mrf,Xfctan'^/3 J='^ 



w 



iw 
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X, 
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Ezl^z. 



2i 72j 



i=i 






■H2{x\ 



1 -^2 )i 



(4.12) 



where x. 



t/c, 



Mf/m, 



C,-5 



1, 2, j = 1,2 are the ratios of the stop and chargino masses 



squared, the matrices Z+ and Z_ are defined in ref. [Q and the function H2{x,y) is 
defined in the Appendix. The factor 

sin^/3 



^z. 



cos^ a 
-Ml 



+ 



sm^ a 



Ml 



T 



M\ 



(4.13) 



depends on the masses of the CP-even neutral Higgs bosons h^ and if°, the mass of the 
CP-odd Higgs boson A^ (in the MSSM M%+ = M\ + M^) and the mixing angles a and 

! 1, sin^ a ^ and Mh ~ Ma- 



(3. For tan/3 > 1 and M^ > 130 GeV, cos^ a 
idR)j {dL)j idL)j 




h'>,H^,A^ 




h'',H'',A^ 




h'^^H^A^ 



Figure 11: Double penguin diagram contributing to: a) Cf"^, b) Cf^^ and c) C^^ Wilson 
coefficients in the MSSM with large tan/3. 



As has been observed in |^, 5^'^'>Cf^^ and (5(o)(7f^^ depend on JF_ which for tan/5 ^ 1 
is close to zero and strongly suppresses these corrections. However the correction ^^^^C^^ 
is proportional to JF,_ which is not suppressed in this limit. Approximating for simplicity 
the dimensionless factor Xtc by unity, it is easy to see that in the case of the B^-Bg mixing 
this correction, although proportional to small strange quark mass {ms{fit) ~ 61 MeV), 
can be for tan/5 ~ 50 and Mh+ ~ 200 GeV as large as (J'-^-'C^^ ~ 2.5 i.e. of the same order 
of magnitude as the SM contribution to C"^^^. Consequently there can be a significant, 
growing as tan^/3, contribution to the C^^ Wilson coefficient which is further enhanced 
(relative to the CY^^ coefficient) by QCD effects (see sec. 3). 
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Figure 12: 1 -|- /^ in the MSSM as a function of the mixing angle of the top squarks for 
different hghter chargino masses and compositions (r = M2/ ^). Sohd, dashed, dotted 
and dot-dashed hues correspond to stop masses (in GeV) (500,650), (500,850), (700,850) 
and (700,1000), respectively. 

An important feature of the double penguin contribution is its fixed negative sign 
(because it is proportional X^c*) i-^- ^^^ same as the sign of the dominant effects of the 
charged Higgs box diagrams at large tan/3. Therefore the double penguin contribution 
interferes destructively with the SM contribution and leads io 1 + fs < 1. Another 
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Figure 13: 1 + /^ in the MSSM for lighter chargino mass 750 GeV, r = M2I [i = —0.5 
and stop masses (in GeV) (500,850), (700,1000), (500,850) and (600,1100) (sohd, dashed, 
dotted and dot-dashed hnes, respectively) as a function of a) tan /3 and b) Mh+ ■ In panel 
a) solid and dashed (dotted and dot-dashed) lines correspond to Mh+ = 200 (600) GeV, 
and in panel a) solid and dashed (dotted and dot-dashed) lines correspond to tan/? = 50 
(35). 



interesting feature is its strong dependence on the left-right mixing of the top squarks 
which is clearly visible in fig. ^ where we show 1 + /^ as a function of the stop mixing 
angle dt for different chargino masses and compositions and different choices of the stop 
masses. For the same value of the mixing angle 6t, larger effects are obtained for bigger 
stop mass splitting because in this case the parameter \At\ has to be larger. It should be 
also stressed that this contribution does not vanish when the mass scale of the sparticles 
is increased (i.e. when all mass parameters are scaled uniformly). Thus, large effects 
decreasing 1 + fs below unity can be present in the MSSM also for the heavy sparticles 
provided the mass scale of the MSSM Higgs sector remains low and tan/5 is large. This 
is illustrated in figs. |l3|a and b where we show 1 + /^ as a function of tan/3 (panel a) 



and Mh+ (panel b). Positive contribution to 1 + /^ seen in fig. |T3|a for tan/3 < 2.5 and 
Mh+ = 200 GeV is due to the ordinary charged Higgs box diagrams which contribute to 
the universal part of fs, fd and fs through the Wilson coefficient of the standard QY^^ 
operator. For lighter H~^ and light charginos 1 + fs can reach values ~ 2 |2^ . As follows 
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from ( |2.18| ) such high values of 1 + /^ could be soon excluded by the measurement of 
AM,. 

The implications of these results are as follows. A large portion of the MSSM param- 
eter space considered here leads to 1 + /, values violating already at present the bound 
( p. IS ) and is, hence, excluded. It is particularly interesting that part of this parameter 



space is not yet excluded by the recent results for the 6 — > S7 rate. This is because the 
theoretical prediction for the latter depends, besides stop and chargino sector parameters 
(on which the double penguin contribution does depend), also on masses of sfermions from 
the first two generations, the gluino mass etc. This can have important consequences for 
other processes involving the B^ mesons. For instance, orders of magnitude enhancement 
of the B^^ ^ fi^fi^ branching ratio found in the MSSM for very large mixing of stops 



55| , pq , |57| , pql is ruled out, at least for sparticles so heavy that their direct contribution 
to the B^-B^ mixing cannot significantly reduce the negative double penguin contribution 
to 1 + fs- Thus, finding BR{B^^ —>■ yU^yU^) close to the present experimental bound and 
no light stops and charginos would strongly indicate a non-minimal flavour violation in 
supersymmetry [Q. 



Furthermore, for the MSSM parameters, for which the bound ( |2.18|) is respected 



1 + /d ~ 1 + /e ~ 1 and in this sector the MSSM in the limit Mw ^ Mh+ < M^particie 
mimics the SM. In particular, it is still consistent with the present experimental data for e, 
AM^, AMg and a^^s- Measuring AM, larger than the present lower bound AMg > 15 /ps 
will further limit allowed combinations of stop mixings, their mass splittings and chargino 
parameters. Only very large values of AM,, requiring 1 + /s ~ 1, would rule out the 
supersymmetric scenario considered here entirely; the SM would be then ruled out too. 

Obviously, finding the asymmetry a^Ks below its SM value would also rule out this 
scenario entirely because for 1 + fd ^ 1 + fe ^ 1 the unitarity of the CKM matrix 
requires sin 2/3 > 0.5, i.e. bigger than the bound ( |1.6|) valid in MVF models which admit 
1 + /d ~ 1 + /e 7^ 1- If CL^Ks is found around 0.6 then the combination of constraints 
from Rf,, e and sin 2(3 similar to the ones shown in fig. |^ can put slightly stronger limits 
on 1 + fs ^ Rsd than the bound ( p.l8| ) alone but the usefulness of this limits will depend 
crucially on the measured value of AM,. 

Finally, if experimentally a^Xs ~ 0.7 and AM, combined with improved lattice results 



for \l Bb^Fb^ and ^ allow for Rsd ~ 1 -|- /, ~ 0.65 — 0.8, this scenario can lead to angle 7 
moderately smaller than the one predicted in the SM. 
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5 Conclusions 

In this paper we have investigated the role of new dimension six four-fermion |AF| = 2 
operators in models with minimal flavour violation. Short distance contributions to the 
mass differences AM^, AM^ and to the CP violation parameter e are parameterized by 
three real functions F^*, F^^ and F^^, respectively. General formulae for F/^ in terms of 
the Wilson coefficients evaluated at the scale /i = /xnp, the relevant QCD renormalization 
group factors and the non-perturbative i?j-factors have been presented. 

We have proposed a few simple strategies involving the ratio AMg/AMa, sin 2/5 and 
the angle 7 that allow to search for the effects of the new operators. We have also 
found model independent bounds on the functions 1 + /j = F/^/S'o(xt) that should be 
considerably improved once AMs/AM^, sin 2/5 and the angle 7 are precisely measured 
and our knowledge about non-perturbative parameters and the CKM elements \Vub\ and 
\Vcb\ is improved. Our findings can be summarized as follows: 

• The present experimental and theoretical uncertainties allow for sizable contribu- 
tions of new operators to AM^^a and e. 

• As the unitarity of the CKM matrix implies \Vts\ ~ \Vcb\ independently of new 
physics contributions, the function 1 + fs can be determined from the experimental 



value of AMg subject to the uncertainties in |Vcfe|, rrit and in particular \J BbsFbs- 
For instance for AMg = 18/ps we find 0.63 < 1 + /s < 1.55. The decrease of the 



theoretical error in \J Bb^Fb^ accompanied by a precise measurement of AM^ should 
tell us whether 1 + /^ > 1 or 1 + /^ < 1 thereby excluding certain scenarios and 
putting important constraints on the parameter space of the surviving models. 

We find that values of Rsd = (1 + /«)/(! + fd) substantially different from unity 
would allow sin 2/3 to be lower than in the MFV models and in particular in the 
SM. Simultaneously Rde = (1 + /d)/(l + /e) < 1 would be favoured. 

Whether Rgd > 1 or Rs^ < 1 is favoured by the data can be decided by the mea- 
surement of the angle 7 with 7 > 90° and 7 < 90° corresponding for AMg = 15/ps 
to Rsd > 1-2 and Rsd < 1-2, respectively. For a given Rsd and sin 2/5, the predicted 
angle 7 decreases with increasing AM^. For Rsd > 1.5 and AMs/AMd < 40 values 
of the angle 7 greater than 90° are possible. 

We have determined the presently allowed ranges in the (1 + Z^, 1 + /g) and {Rsd, 1 + 
fs) planes. An analysis of a hypothetical measurements of AMs and a^K^ that allow 
the determination of the unitarity triangle illustrated various possibilities further. 
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As an example we have analyzed the role of new operators in the MSSM with large 
tan/3 = f2/fi in the limit of heavy sparticles, investigating in particular the impact of 
the extended Higgs sector on the unitarity triangle. Here our findings are as follows: 

• The largest effects of new contributions for large tan/3 are seen in AM^. The 
corresponding contributions to AM^; and e are strongly suppressed either by inverse 
powers of tan (3 or by the smallness of d-quark mass. 

• The dominant contributions to AM^ for large tan (3 come from the operator Q\^ = 
(6(1 — 75)5) (6(1 + 75)3). They originate from double penguin diagrams involving 
neutral Higgs particles and, to a lesser extent, in the box diagrams with charged 
Higgs exchanges. The dominant double penguin diagrams arise through the gener- 
ation of flavour non-diagonal tan/3 enhanced couplings of neutral Higgs bosons to 
the down-type quarks and depend strongly on the mixing of the top squarks and 
their mass splitting. 

• The contribution of double penguins grows like tan'^ /3 and interferes destructively 
with the SM contribution, suppressing considerably I + fs below unity. 

• All these findings have the following phenomenological consequences. The MSSM 
with large tan/?, substantial stop mixing and large stop mass splitting realizes the 
Rsd < 1 scenario with 7 < 90° and generally smaller than in the SM and MFV 
models. As Rde = 1 and 1 + /^ = 1, the lower bound sin 2/3 > 0.50 valid in the 
SM remains unchanged. Consequently if a^Xs is found below 0.50, this scenario of 
the MSSM will be excluded (together with the SM) while other MSSM scenarios 
with lighter sparticles and lower tan/3, belonging to the MFV class, may still be 
consistent with the data. As seen in fig. ^, for higher values of a^Ks the MSSM 
scenario considered here is a vital possibility with the angle 7 smaller than in the 
SM, although values of Rsd as low as 0.6 appear rather improbable. 



• The constraint ( |2.18| ), which basically limits the magnitude of the stop mixing pa- 
rameter At, has also important consequences for other processes involving the B^ 
mesons. For example, in the scenario considered in this paper, it severely limits 
possible enhancement of the B^^ ^ fJ'^fJ'^ decay rate. 

Detailed analysis of AMd^s and e in the MSSM at large tan/3, including also scenarios 
with light sparticles will be presented soon ^6 . 



It will be exciting to watch future developments in the experimental values of AM^, 
a^Ks ^iid the angle 7 that will either choose one of the possibilities considered in this 
paper or constrain the parameters of GMFV models. 
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Appendix A 

The four point loop functions with zero external momenta are defined as follows: 

n( h ri\ = /" ^^^ ^ ^ 1 ^ 

Uo[a,o,c,a) J ^2 ^a][b][c][d\ (h - a){h - c){h - d)^^ a 

1 l/'A-l\ 

(c ~ a){c — h) (c — d) a {d — a) {d — b){d — c) a 

Hi h ri\ - M'^ '^' ^' 1 ^ 

1^2 la, 0, c, a J j ^2 ^ [5] [c] [rf] {h-a){h-c){h-d) ^^ a 

+ 7 77 TTT 3Tlog- + 73 T73 ITTl T ^^S " (^-2) 

(c — a j (c — 0) [c — d) a [d — a) [d — b)[d ~ c) a 

where [a] = k'^ — a etc. Finally, 

TT / \ x\nx ylny / A o\ 

H2{x, y) = r + . A.3 

{l-x){x~y) {l-y){y-x) 
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